Rotating solenoidal perfect fluids of Petrov type D. 
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Abstract. We prove that aligned Petrov type D perfect fluids for which the vorticity 
O ' vector is not orthogonal to the plane of repeated principal null directions and for which 

| /' | ■ the magnetic part of the Weyl tensor with respect to the fluid velocity has vanishing 

\ divergence, are necessarily purely electric or locally rotationally symmetric. The LRS 

metrics are presented explicitly. 
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1. Introduction 

Recently there has been some interest in perfect fluid solutions of the Einstein field 
equations, 

i 

G ab = Rab ~ -RQab = (w + p)u a U b + V9ab, (1) 

which have been called [H[2] solenoidal. These are space-times having a vanishing spatial 
divergence of the magnetic part of the Weyl tensor, when the latter is decomposed [3] 
with respect to the fluid velocity u: 

(divH) a = D b H ab = h a m h bn H mb . n = 0. (2) 

The resulting models form an obvious generalization of the purely electric (H = 0) 
fluids, among which one finds many physically relevant exact solutions. They also offer 
the hope of being mathematically tractable, as one set of the Bianchi equations pi] 
simplifies considerably and becomes algebraic in E, <t, uj and w + p: 

3E ab uj b +(w+ p)u a - [<r, E} a = 0. (3) 

All investigations so far seem to indicate however that the solenoidal condition - 
although significantly weaker than the condition H = — is still very strong and leads 
to classes of well known solutions, such as the locally rotationally symmetric (LRS) ones 
or the purely electric solutions. In the present paper we show that if an aligned perfect 
fluid of Petrov type D (meaning that the fluid velocity is in the plane S of the repeated 
principal null directions k and 1), is solenoidal and has a vorticity vector which is not 
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orthogonal to E, then this fluid is necessarily purely electric or LRS. In the latter case 
the line-elements can be given in explicit form. 

We will make use of the Newman-Penrose formalism and follow the notation and 
sign conventions of |3j, whereby the Newman-Penrose equations (7.21a - 7.21r) and 
Bianchi identities (7.32a - 7.32k) will be indicated as (npl - npl8) and (bl - bll) 
respectively. We choose a canonical type D tetrad, in which = <fi + iip is the only 
non- vanishing component of the Weyl-spinor (0, ip real). Provided that the fluid velocity 
is aligned with E, we fix a boost such that the only non- vanishing components of the 
traceless Ricci spinor are given by $ o — ^22 = 2$n = 5/4, with S = w +p the sum of 
the fluid's energy density and pressure. The energy density is then related to the Ricci 
scalar by the relation R = Aw — 3S. We also decompose the fluid's vorticity vector u) 
in a component uj± orthogonal to E and a component parallel to E. One can easily 
verify, see for example [5 J , that u) — uj±_ is proportional to the following combination of 
the spin-coefficients, p — p + ~p — p, which we will henceforth abbreviate as 

uj^ = p-p + p-p (4) 

and which we assume to be non-zero. 



2. Main result 

Choosing a canonical type D tetrad as described above, the solenoidal condition (j2J) can 

be expressed by the equations 
3 

6lp = ~(V — 7T + T — K)ljj (5) 

(D-A)^ = |(p + p + /i + 7l)V. (6) 
The Bianchi equations on the other hand lead to the differential equations 

5S = -12K((j) + t^) + (K-W + 2f3 + 2a)S (7) 

1 i 
5(f) = -(2k + 3r - 3tT)0 + -(2k + 3r + 37f)V> 

S 

- — (K + V-T-W + Afi + Aa) (8) 

3 3 S 

dip = —-i(r + 7f)0 + -(t — W)if) - i— (k + V - t - W) (9) 
2 2 8 

(D - Aty = \{p + p + P + 7i# + {MS + 60) (10) 

3 3 
(D + A)ip = -(p + p - p -p)ip + -i(p - p + p -p)(j) (11) 

and three algebraic relations, obtained by making linear combinations of (b2, b5, b7,b8), 
namely 

S 

3(k + V)<f) + 3z(k - z7)-0 - -(k + V - t - 7f + 4a + A/3) = (12) 

12a{(j> + + S{a - A) = (13) 
a{<f) + iip) + A(0-z'V>) = (14) 
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Comparing flS]) and (fTDl) we see that u^(S + 60) = and hence, provided that u>z ^ 0, 

5 + 60 = 0. (15) 

Herewith the determinant of the linear and homogeneous system (fl3l [141) in a and A 
equals Gtjj 2 . Discarding for the moment!]]] the purely electric case (as then the condition 
divH = is trivially satisfied), it follows that the null congruences corresponding to k 
and 1 are shear-free: 

ct = A = 0. (16) 

In addition the equations (jEJ), (EJ) and ( TT2l) imply 

/3 + a = r — K + n — V (17) 

and 

(k + V - 3r - 3vf)0 - 2x(k - 17)^ = 0. (18) 

A second algebraic relation between k, r, 17 and 7f is obtained by applying the 5-operator 
to ([15]): using (J7JH]) this yields 

(3k - 517 - r + 77f)0 + 2i(2n - 3r - 3tt)^ = 0. (19) 
Solving ffl8l) and ffl9|) for k and r one finds 

k = [-60 2 vf + (2i(pip + 40 2 - 3^ 2 )^] (2iH + 20 2 - 3^ 2 ) ~* (20) 

r = [(2i(j)ip - 40 2 + 3^ 2 )vf + 2(0 2 - ip 2 )v] {2i<pi) + 20 2 - 3^ 2 ) ^ . (21) 

Next we notice that the NP-equations np2, npl6 and the complex conjugates of np7, 
nplO give algebraic expressions for the (^-derivatives of k,t,u and IT. Substituting in 
these the two expressions (1201 |2T]) . two further restrictions on 7r and v follow: 

3(3^ 2 - 20 2 - 10^0) (0 - ^)0V + 0(80 4 + 28i0 3 ^ + 38^ 2 2 - 21^ 4 - 32z0^ 3 )z/tt 
-(20 5 + 12i0> - 23^ 4 + 12# 5 + 210V - 33i0 2 ^ 3 )z/ 2 = 

0(80 5 + 84i0 4 ^ + 1020V + 54^ 5 + 81^ 4 + 44i0 2 ^ 3 )vr 2 + 6i(-160 5 + 26i0 4 ^ 
-12^ 4 + 180 3 ^ 2 - Qiip 5 - 5i0V)^vr + (94^ 2 4 + 12^ 6 - 80 6 
-112V0 3 + 44#0 5 + 66# 5 - 99^ 4 2 )z/ 2 = 0. 

Elimination of v results then in 

(1440 s + 424^ 2 6 + 848^ 4 4 + 838^ 6 2 + 297^ 8 )(0 + 2iip) 2 
x(3ip 2 - 20 2 + 2^0) 4 2 V> 2 7T 4 = 0, 

while elimination of it gives a similar equation for v. As the first factor is positive 
definite, it follows that solutions are purely electric or purely magnetic, with the second 
option being excluded because of ([3]), or satisfy v = n = 0. In the latter case n = r = 
follow from (12011211) and solutions are LRS (see [5]), more particularly of Ellis' class I 

mm- 

II a detailed discussion of the purely electric Petrov type D perfect fluids will be given in a forthcoming 
publication, see [6] 
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3. Line elements and discussion 

In this section we present the line elements of the LRS I aligned perfect fluid space-times 
of Petrov type D satisfying the solenoidal condition divH = and having vorticity not 
orthogonal to E. As the general form of the rotating LRS class I perfect fluid line 
elements is well known [7j, it suffices to use an algebraic computer package such as 
GRTensor [9J in order to explicitly express the solenoidal condition. It turns out however 
that the following metric form is more convenient for the actual computations: 

ds * = _0 O 2 + 012 + ^2 + ^ (22) 

with the duals l of the basisvectors ej given by 
0° = f~ l s{dt + qr 2 {l + ^-Y l d V ] 

e l = r 2 dx 

e 2 = f-\i + ^dr 

3 = rf- 1 (l + —)- 1 d V (23) 

(/, s functions of x and k, q constants). This represents a perfect fluid with velocity e 
provided the equation of pressure isotropy holds, 

f 2 s" - 2s ff" + fs'f - 2sf' 2 + ks + 2q 2 s 3 = 0. (24) 

The eo congruence is non-expanding and shearfree, with vorticity and acceleration 
respectively given by 

u) = qfse±, ii = -sf(f/s)'e l . (25) 

The solutions therefore belong to Collins' class IIIAGii of stationary and shear-free 
fluids obeying a barotropic equation of state and having vorticity parallel to the 
acceleration [10]. They are of class LRS Id in the classification of Stewart and Ellis [8] 
and admit a isometry group acting multiply transitively on the timelike hypersurfaces 
orthogonal to e±. Calculating the electric and magnetic part of the Weyl tensor results 
in 



H = g/Vdiag(0,2,-1,-1), (26) 

P 
6s 



E = f {f 2 s" + ffs' - Aq 2 s 3 - ks) diag(0, 2, -1, -1). (27) 



while 



divH = 2qfs"e l . (28) 

The solenoidal condition therefore holds when s = 1 or s = x. In the first case, s = 1, 
solutions are purely electric generalizations of the Godel metric, with / given by 

/= ( Cl + C2X+ (^ + g 2 )x 2 ) l/2 ; (2Q) 
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while in the second case, s — x, the solution / of (1211) is given by 

/ = ( Cl + C2X 3/2 + kx 2 + ^4)1/2 (30) 

(ci, c 2 constants of integration). In the latter case the expressions for the magnetic and 
electric parts of the Weyl curvature reduce to 

H = g/ 3 diag(0,2,-l,-l), (31) 
E = ^-(5ci - 5c 2 x 3/2 + q 2 x 4 + 5fcr 2 )(24gV + 5c 2 x 1/2 ) diag(0, 2, -1, -1). (32) 

Two simple examples are obtained for c\ = c 2 = and c\ = 0, c 2 = in which the 
pressure p and matter density w are related respectively by the expressions 

1 18 2 4 3 4 6 . . 

p = -w kq x = — q x (33) 

y 5 25 y 25 y V 7 

and 

P = - ^(5A: + 2g 2 x 2 ) 2 = ^(gV ~ 5k)(2q 2 x 2 + 5k) 2 . (34) 

In the flat case, k — 0, both metrics reduce to 

ds 2 = 5g~ 2 aT 2 (-(dt + qr 2 d<p) 2 + 5<r 2 aT 6 cfe 2 + x" 2 (rfr 2 + r 2 d0 2 )) . (35) 
Introducing standard coordinates by putting x = y/E/ (qz) this can also be written as 

2 4 

ds 2 = -z 2 (dt + qr 2 dip) 2 + — ds£, (36) 



with ds 2 = dr 2 + r 2 d(p 2 + <iz 2 the metric of 3-d Euclidean space. The solution clearly 
satisfies then the regularity condition on the r = axis and represents a rigidly rotating 
perfect fluid obeying a 7-law equation of state p = ~w. This particular equation of 
state suggests a possible link with the Collins-Stewart [11] purely magnetic (spatially 
homogeneous and expanding, but non-rotating) perfect fluid metric 

ds l = _ a 2 d f2 + f/^ dx + r 2 rf0) 2 + t 4/3 (rfr 2 + ^2^ ^2 = 1 (37) 

o 

However, although a complex coordinate transformation t — > gx, x — > iq^ 1 ^ 1 ^ 3 
transforms (I3"5]l into a metric of the form (1371) (modulo a global scale factor), the resulting 
constant a 2 does not have the correct value a 2 = |. The metric f l36|) belongs to a class 
of self-similar solutions with an H 5 similarity group discussed by Wainwright in [T2] . 
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